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Colligative properties of solutions: 

On . 

o II. Vanishing concentrations 

(N 

^ ■ Kenneth S. Alexander,^ Marek Biskup,^ and Lincoln Chayes^ 



We continue our study of colligative properties of solutions initiated in ref. [1]. 
I We focus on the situations where, in a system of linear size L, the concentration 

• and the chemical potential scale like c = ^/L and h = b/L, respectively. We 

I find that there exists a critical value such that no phase separation occurs 

■ for £, < £,t while, for ^ > ^t, the two phases of the solvent coexist for an interval 
I of values of b. Moreover, phase separation begins abruptly in the sense that a 

macroscopic fraction of the system suddenly freezes (or melts) forming a crystal 
(or droplet) of the complementary phase when b reaches a critical value. For 

■ certain values of system parameters, under "frozen" boundary conditions, phase 

separation also ends abruptly in the sense that the equilibrium droplet grows 

■ continuously with increasing b and then suddenly jumps in size to subsume the 
I entire system. Our findings indicate that the onset of freezing-point depression 

is in fact a surface phenomenon. 
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1. INTRODUCTION 
1.1. Overview 



In a previous paper (ref. [1], henceforth referred to as Part I) we defined a model 
of non-volatile solutions and studied its behavior under the conditions when 
the solvent undergoes a liquid-solid phase transition. A particular example of 
interest is the solution of salt in water at temperatures near the freezing point. 
In accord with Part I we will refer to the solute as salt and to the two phases of 
solvent as ice and liquid water. 

After some reformulation the model is reduced to the Ising model coupled 
to an extra collection of variables representing the salt. The (formal) Hamilto- 
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nian is given by 
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Here we are confined to the sites of the hypercubic lattice Z'^ with d > 2, 
the variable ax e {+1, —1} marks the presence of liquid water (ax — 1) 
and ice (ax — —1) at site x, while G {0, 1} distinguishes whether salt is 
present (S^ = 1) or absent (S^ = 0) at x. The coupling between the a's is 
ferromagnetic (J > 0), the coupling between the a's and the S's favors salt in 
liquid water, i.e., k > — this reflects the fact that there is an energetic penalty 
for salt inserted into the crystal structure of ice. 

A statistical ensemble of direct physical — and mathematical — relevance 
is that with fluctuating magnetization (grand canonical spin variables) and a 
fixed amount of salt (canonical salt variables). The principal parameters of the 
system are thus the salt concentration c and the external field h. As was shown 
in Part I for this setup, there is a non-trivial region in the (c, /i) -plane where 
phase separation occurs on a macroscopic scale. Specifically, for (c, h) in this 
region, a droplet which takes a non-trivial (i.e., non-zero and non-one) fraction 
of the entire volume appears in the system. (For "liquid" boundary conditions, 
the droplet is actually an ice crystal.) In "magnetic" terms, for each h there 
is a unique value of the magnetization which is achieved by keeping part of 
the system in the liquid, i.e., the plus Ising state, and part in the solid, i.e., the 
minus Ising state. This is in sharp contrast to what happens in the unperturbed 
Ising model where a single value of h (namely, h = 0) corresponds to a whole 
interval of possible magnetizations. 

The main objective of the present paper is to investigate the limit of in- 
finitesimal salt concentrations. We will take this to mean the following: In 
a system of linear size L we will consider the above "mixed" ensemble with 
concentration c and external field h scaling to zero as the size of the system, L, 
tends to infinity. The goal is to describe the asymptotic properties of the typical 
spin configurations, particularly with regards to the formation of droplets. The 
salt marginal will now be of no interest because salt particles are so sparse that 
any local observable will eventually report that there is no salt at all. 

The main conclusions of this work are summarized as follows. First, in 
a regular system of volume V = L'^ oi characteristic dimension L, the scal- 
ing for both the salt concentration and external field is L^. In particular, we 
should write h — bL~^ and c — iL~^. Second, considering such a system 
with boundary condition favoring the liquid state and with h and c enjoying the 
abovementioned scalings, one of three things will happen as ^ sweeps from 
to infinity: 
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(1) If 6 is sufficiently small negative, the system is always in the liquid state. 

(2) If b is of intermediate (negative) values, there is a transition, at some {(6) 
from the ice state to the liquid state. 

(3) Most dramatically, for larger (negative) values of b, there is a region — 
parametrized by < C < ^2{b) — where (macroscopic) phase sepa- 
ration occurs. Specifically, the system holds a large crystalline chunk of 
ice, whose volume fraction varies from unity to some positive amount 
as ^ varies from ^i{b) to ^2{b). At ^ = ^2{b), all of the remaining ice 
suddenly melts. 

We obtain analogous results when the boundary conditions favor the ice state, 
with the ice crystal replaced by a liquid "brine pocket." However, here a new 
phenomenon occurs: For certain choices of system parameters, the (growing) 
volume fraction occupied by the brine pocket remains bounded away from one 
as ^ increases from ^i{b) to ^2{b), and then jumps discontinuously to one at 
^2{b). In particular, there are two droplet transitions, see Fig. 1. 

Thus, we claim that the onset of freezing point depression is, in fact, a 
surface phenomenon. Indeed, for very weak solutions, the bulk behavior of 
the system is determined by a delicate balance between surface order devia- 
tions of the temperature and salt concentrations. In somewhat poetic terms, the 
predictions of this work are that at the liquid-ice coexistence temperature it is 
possible to melt a substantial portion of the ice via a pinch of salt whose size 
is only of the order V^~d. (However, we make no claims as to how long one 
would have to wait in order to observe this phenomenon.) 

The remainder of this paper is organized as follows. In the next section 
we reiterate the basic setup of our model and introduce some further objects of 
relevance. The main results are stated in Sections 2.1-2.3; the corresponding 
proofs come in Section 3. In order to keep the section and formula numbering 
independent of Part I; we will prefix the numbers from Part I by "1." 

1 .2. Basic objects 

We begin by a quick reminder of the model; further details and motivation are 
to be found in Part 1. Let A C Z"' be a finite set and let dA denote its (external) 
boundary. For each x E A, we introduce the water and salt variables, G 
{ — 1, +1} and Sx G {0, 1}; on 9A we will consider a fixed configuration aoA G 
{— 1, +1}^^. The finite-volume Hamiltonian is then a function of (cta, Sa) and 
the boundary condition agA that takes the form 

pJ^A{aA,SA\aoA) = - J ^ a^ay - h'^ax + k'^Sx — (1.2) 

{x,y) xeA xeA 

xeA,yeZ'' 
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Fig. 1. The phase diagram of the ice-water system with Hamihonian (1.1) and fixed salt 
concentration c in a Wulff-shaped vessel of linear size L. The left plot corresponds to the 
system with plus boundary conditions, concentration c = ^/L and field parameter h = b/L, 
the plot on the right depicts the situation for minus boundary conditions. It is noted that as ^ 
ranges in (0, oo) with b fixed, three distinct modes of behavior emerge, in the L ^ cx) limit, 
depending on the value of b. The thick black lines mark the phase boundaries where a droplet 
transition occurs; on the white lines the fraction of liquid (or solid) in the system changes 
continuously. 



Here, as usual, {x, y) denotes a nearest-neighbor pair on 7/ and the parame- 
ters J, K and h represent the chemical affinity of water to water, negative affin- 
ity of salt to ice and the difference of the chemical potentials for liquid-water 
and ice, respectively. 

The a priori probability distribution of the pair (cta, Sa) takes the usual 
Gibbs-Boltzmann form PX^''((Ta, Sa) oc e-^-^^^^^'^^l^s^^. For reasons ex- 
plained in Part I, we will focus our attention on the ensemble with a fixed 
total amount of salt. The relevant quantity is defined by 

Nk = Y,S.. (1.3) 
The main object of interest in this paper is then the conditional measure 

^r'^'X-) = i^A = Lc|Aij), (1.4) 

where |A| denotes the number of sites in A. We will mostly focus on the situ- 
ations when (T9A = +1 or asK = —1, i.e., the plus or minus boundary condi- 
tions. In these cases we denote the above measure by P^''^'^, respectively. 
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The surface nature of the macroscopic phase separation — namely, the 
cases when the concentration scales Uke the inverse linear scale of the system — 
indicates that the quantitative aspects of the analysis may depend sensitively on 
the shape of the volume in which the model is studied. Thus, to keep this work 
manageable, we will restrict our rigorous treatment of these cases to volumes 
of a particular shape in which the droplet cost is the same as in infinite volume. 
The obvious advantage of this restriction is the possibility of explicit calcula- 
tions; the disadvantage is that the shape actually depends on the value of the 
coupling constant J. Notwithstanding, we expect that all of our findings are 
qualitatively correct even in rectangular volumes but that cannot be guaranteed 
without a fair amount of extra work; see [17] for an example. 

Let V C M*^ be a connected set with connected complement and unit 
Lebesgue volume. We will consider a sequence (Vl) of lattice volumes which 
are just discretized blow-ups of V by scale factor L: 

VL^{xeZ'^:x/LeV}. (1.5) 

(The sequence of L x • • • x L boxes (A^^) from Part I is recovered by let- 
ting V — [0, 1)^.) The particular "shape" V for which we will prove the 
macroscopic phase separation coincides with that of an equilibrium droplet — 
the Wulff-shaped volume — which we will define next. We will stay rather suc- 
cinct; details and proofs can be found in standard literature on Wulff construc- 
tion ( [2,4, 5, 7, 8, 12] or the review [6]). Readers familiar with these concepts 
may consider skipping the rest of this section and passing directly to the state- 
ments of our main results. 

Consider the ferromagnetic Ising model at coupling J > and zero ex- 
ternal field and let P^''^ denote the corresponding Gibbs measure in finite vol- 
ume A C Z'^ and plus/minus boundary conditions. As is well known, there 
exists a number — Jc{d), with Jc(l) = oo and Jc{d) e (0, oo) if d > 2, 
such that for every J > the expectation of any spin in A with respect to Pj''^ 
is bounded away from zero uniformly in A C Z'^. The limiting value of this ex- 
pectation in the plus state — typically called the spontaneous magnetization — 
will be denoted by = mi,{J). (Note that m^^ = for J < J^. while > 
for J > Jc.) 

Next we will recall the basic setup for the analysis of surface phenomena. 
For each unit vector n E W^, we first define the surface free energy rj(n) 
in direction n. To this end let us consider a rectangular box V{N, M) C 
with "square" base of side N and height M oriented such that n is orthogonal 
to the base. The box is centered at the origin. We let Z^'^^ denote the Ising 
partition function in V{N, M) fl Z'^ with plus boundary conditions. We will 
also consider the inclined Dobrushin boundary condition which takes value +1 
at the sites x of the boundary oiV{N, M) n for which x ■ n > {) and — 1 
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at the other sites. Denoting the corresponding partition function by Zpj^, the 
surface free energy Tj(n) is then defined by 




The Umit exists by subadditivity arguments. The quantity Tj(n) determines 
the cost of an interface orthogonal to vector n. 

As expected, as soon as J > Jc, the function n i— Tj{n) is uniformly 
positive [14]. In order to evaluate the cost of a curved interface, rj(n) will 
have to be integrated over the surface. Explicitly, we will let J > Jc and, given 
a bounded set y c R*^ with piecewise smooth boundary, we define the Wulff 
functional Wj by the integral 



where AA is the (Hausdorff) surface measure and n is the position-dependent 
unit normal vector to the surface. The Wuljf shape W is the unique minimizer 
(modulo translation) ofV ^ W{V) among bounded sets F c M'' with piece- 
wise smooth boundary and unit Lebesgue volume. We let {Wl) denote the 
sequence of Wulff-shaped lattice volumes defined from V — W via (1.5). 



We are now in a position to state and prove our main results. As indicated 
before, we will focus on the limit of infinitesimal concentrations (and external 
fields) where c and h scale as the reciprocal of the linear size of the system. 
Our results come in four theorems: In Theorem 2.1 we state the basic surface- 
order large-deviation principle. Theorems 2.2 and 2.3 describe the minimizers 
of the requisite rate functions for liquid and ice boundary conditions, respec- 
tively. Finally, Theorem 2.4 provides some control of the spin marginal of the 
corresponding Gibbs measure. 

2.1. Large deviation principle for magnetization 

The control of the regime under consideration involves the surface-order large- 
deviation principle for the total magnetization in the Ising model. In a finite 
set A C Z'^, the quantity under considerations is given by 




(1.7) 



2. MAIN RESULTS 




(2.1) 
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Unfortunately, the rigorous results available at present for d > 3 do not cover 
all of the cases to which our analysis might apply. In order to reduce the amount 
of necessary provisos in the statement of the theorems, we will formulate the 
relevant properties as an assumption: 

Assumption A Let d > 2 and let us consider a sequence of Wulff-shape 
volumes Wl- Let J > Jc and recall that denotes the Gibbs state of the 
Ising model in Wl, with ±-boundary condition and coupling constant J. Let 
rrii, = 'm^,{J) denote the spontaneous magnetization. Then there exist functions 
■^±,j '■ [—^-k-, rrii] — >• [0, oo) such that 

Um lim -J— logP^'-^(|Mi - mL'^l < eL'^) = Am) (2.2) 

holds for each m e [— rUi]. Moreover, there is a constant wi e (0, oo) 
such that 

^±,j(m) = [-^^ ) ^1 (2-3) 
is true for all m e [— m^, m*] . 

The first part of Assumption A — the surface-order large-deviation princi- 
ple (2.2) — has rigorously been verified for square boxes (and magnetizations 
near ±m^) in d — 2 [8, 12] and in d > 3 [5, 7]. The extension to Wulff- 
shape domains for all m e [— my,,myf] requires only minor modifications in 
d = 2 [16]. For d > 3 Wulff-shape domains should be analogously control- 
lable but explicit details have not appeared. The fact (proved in [16] for d = 2) 
that the rate function is given by (2.3) for all magnetizations in [— m*, m*] is 
specific to the Wulff-shape domains; for other domains one expects the formula 
to be true only when jm^, T "7,| is small enough to ensure that the appropriately- 
sized Wulff-shape droplet fits inside the enclosing volume. Thus, Assump- 
tion A is a proven fact for d = 2, and it is imminently provable for d > 3. 

The underlying reason why (2.2) holds is the existence of multiple states. 
Indeed, to achieve the magnetization m e (— m*, m^) one does not have to alter 
the local distribution of the spin configurations (which is what has to be done 
for m ^ [— m^, mj,]); it suffices to create a droplet of one phase inside the other. 
The cost is just the surface free energy of the droplet; the best possible droplet is 
obtained by optimizing the Wulff functional (1.7). This is the content of (2.3). 
However, the droplet is confined to a finite set and, once it becomes sufficiently 
large, the shape of the enclosing volume becomes relevant. In generic volumes 
the presence of this additional constraint in the variational problem actually 
makes the resulting cost larger than (2.3) — which represents the cost of an 
unconstrained droplet. But, in Wulff-shape volumes, (2.3) holds regardless of 
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the droplet size as long as |m| < m^. An explicit formula for ^±,j(m) for 
square volumes has been obtained in d — 2 [17]; the situation in d > 3 has 
been addressed in [10, 11]. 

On the basis of the above assumptions, we are ready to state our first main 
result concerning the measure P^'^'^ with c ^/L and h ~ b/L. Using 9 
to denote the fraction of salt on the plus spins, we begin by introducing the 
relevant entropy function 

T(m, 9) = -9 log - (1 - ^) log (2.4) 

1 + m 1 — m 

We remark that if we write a full expression for the bulk entropy, S(m, 9; c), 
see formula (3.5), at fixed m, c and 9, then, modulo some irrelevant terms, 
the quantity T(m, 9) is given by {d/dc)E{m, 9; c) at c = 0. Thus, when we 
scale c ~ ^/L, the quantity ^T{m,9) represents the relevant (surface order) en- 
tropy of salt with m and 9 fixed. The following is an analogue of Theorem 1.2. 1 
from Part I for the case at hand: 

Theorem 2.1. Let d > 2 and let J > Jc{d) and k > be fixed. Let 

rrii, — m^,{J) denote the spontaneous magnetization of the Ising model. Sup- 
pose that (2.2) in Assumption A holds and let {cl) and {hi) be two sequences 
such that Cl > Ofor all L and that the limits 

^ = lim Lcl and h = lim LHl (2.5) 

L— >oo -L— >oo 

exist and are finite. Then for all m e [— m*,m*], 

lim lim -^logP^;'^'''^(|ML - mL'^l < eL"^) 

6j,0 L/ — >oo j_j 

= -g±(m)+ inf g±^(mO, (2.6) 

where Q^^{m) — iniQ^^^Q ^ ^^^{m, 9) with 

^±^(m, 9) = -hm - ^k9 - ^T{m, 9) + ^±,j(m), (2.7) 

Various calculations in the future will require a somewhat more ex- 
plicit expression for the rate function m Q^^{m) on the right-hand side 
of (2.6). To derive such an expression, we first note that the minimizer 
of 6* H- > =S^^(m, 9) is uniquely determined by the equation 



9 1 + m ^ 

— e 

1-9 1-m 



(2.8) 
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Plugging this into ^^^{m, 9) tells us that 

QtM) = -^"^ - ^9{m) + ^±,j(m), (2.9) 

where 

,(m) = log(i^ + e"l±^). (2.10) 
Clearly, g is strictly concave for any k > 0. 

2.2. Macroscopic phase separation — "liquid" boundary conditions 

While Theorem 1.2. 1 of Part I and Theorem 2. 1 above may appear formally 
similar, the solutions of the associated variational problems are rather different. 
Indeed, unlike the "bulk" rate function Gh,c{'inn) of Part I, the functions Q^^{m) 
are not generically strictly convex which in turns leads to a possibility of having 
more than one minimizing m. We consider first the case of plus (that is, liquid 
water) boundary conditions. 

Let d > 2 and let J > Jc{d) and k > be fixed. To make our formulas 
manageable, for any function : [— m^, m*] — > R let us use the abbreviation 

_ H^^) - 0(-^*) 

for the slope of (f) between —m^ and m*. Further, let us introduce the quantity 

and the piecewise linear function 62 : [0, 00) — > R which is defined by 

h(^)={ T,'^"'" (2-13) 

Our next result is as follows: 

Theorem 2.2. Let d > 2 and let J > Jc{d) and k > be fixed. Let 
the objects Q'^^, ^nd 62 be as defined above. Then there exists a (strictly) 
decreasing and continuous fiinction hi : [0, 00) IR with the following prop- 
erties: 

(1) hiO > hiOforain > 0, andhiO = HO < 6- 

(2) b[ is continuous on [0, 00), b[{^) — > —g'{mi,) as ^ — > 00 and hi is 
strictly convex on [^t; oo)- 
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(3) For b ^ bi{^),b2{C), the function m ^ Q'l^{m) is minimized by a single 
number m — m+(6, ^) e [—m^-i which satisfies 

{=m^, ifb>bi{^), 
e(-m.,m.), ifb2{0<b<bi{0, (2-14) 

= -m*, ifb<b2(0- 

(4) The function b m^{b,^) is strictly increasing for b G [62(0) (01' 
is continuous on the portion of the line b = b2{^) for which C, > and 
has a jump discontinuity along the line defined by b = bi{^). The only 
minimizers at b — bi{^) and b — &2(0 '^^^ the corresponding limits 
ofb ^ m+{b,^). 

The previous statement essentially characterizes the phase diagram for the 
cases described in (2.5). Focusing on the plus boundary condition we have the 
following facts: For reduced concentrations ^ exceeding the critical value ,^t, 
there exists a range of reduced magnetic fields b where a non-trivial droplet ap- 
pears in the system. This range is enclosed by two curves which are the graphs 
of functions foi and 62 above. For b decreasing to bi{^), the system is in the pure 
plus — i.e., liquid — phase but, interestingly, at bi a macroscopic droplet — an ice 
crystal — suddenly appears in the system. As b further decreases the ice crys- 
tal keeps growing to subsume the entire system when b = 62(0- ^ 
no phase separation occurs; the transition at 6 = &i(0 — ^2(0 is directly 
from m = torn — — m^. 

It is noted that the situation for ^ near zero corresponds to the Ising model 
with negative external field proportional to In two-dimensional setting, 
the latter problem has been studied in [16]. As already mentioned, the gener- 
alizations to rectangular boxes will require a non-trivial amount of extra work. 
For the unadorned Ising model (i.e., c — 0) this has been carried out in great 
detail in [17] for d = 2 (see also [13]) and in less detail in general dimen- 
sions [10,11]. 

It is reassuring to observe that the above results mesh favorably with 
the corresponding asymptotic of Part I. For finite concentrations and exter- 
nal fields, there are two curves, c 1— > /i+(c) and c 1-^ h-{c), which mark the 
boundaries of the phase separation region against the liquid and ice regions, 
respectively. The curve c 1— > /i+(c) is given by the equation 

Mc) = ^log[^, (2.15) 

where (g'+, g_) is the (unique) solution of 

Oa^ 1 + 1 — 
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The curve c i— > /i-(c) is defined by the same equations with the roles of 

and —rrii, interchanged. Since h±{0) = 0, these can be linearized around the 
point (0, 0). Specifically, plugging b/L for h and ^/L for c into h — h±{c) and 
letting L — > oo yields the linearized versions 

6± = /i'±(0)e (2.17) 

of h+ and h^. It is easy to check that h'^{0) = — (^'(im^) and so, in the 
limit ^ — s> oo, the linear function 6+ has the same slope as bi while 6_ has the 
same slope as 62 above. Theorem 2.2 gives a detailed description of how these 
linearized curves ought to be continued into (infinitesimal) neighborhoods of 
size 1/L around (0, 0). 

2.3. Macroscopic phase separation— "ice" boundary conditions 

Next we consider minus (ice) boundary conditions, where the requisite liquid 
water, phase separation and ice regions will be defined using the functions 
61 > 62- As for the plus boundary conditions, there is a value > where the 
phase separation region begins, but now we have a new phenomenon: For some 
(but not all) choices of J and k, there exists a nonempty interval (^t, ^u) of ^ 
for which two distinct droplet transitions occur. Specifically, as b increases, the 
volume fraction occupied by the droplet first jumps discontinuously at 62(0 
from zero to a strictly positive value, then increases but stays bounded away 
from one, and then, at 6 = jumps discontinuously to one; i.e., the ice 
surrounding the droplet suddenly melts. 

For each J > Jc{d) and each k,, consider the auxiliary quantities 

(Note that, due to the concavity property of g, both ^1 and ^2 are finite and 
positive.) The following is a precise statement of the above: 

Theorem 2.3. Let d > 2 and let J > Jc{d) and k > be fixed. Then 

there exist two (strictly) decreasing and continuous functions 61, 62 : [0, 00) — > 
R and numbers ^t,Cu G (0, 00) with < ■s'mc/i that the following properties 
hold: 

(1) hiO > hiOforall^y 0, andhiO = HO iffi < ft- 

(2) 62 is strictly concave on [^tj b^i^i) — > —g'{—m^) as ^ — > 00, 61 is 
strictly convex on (^t) Cu) (^nd, outside this interval, 

i.(i)={f~^'^'" (2.19) 
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(3a) If^i > ^2. then = = <3nJ 63 continuous on [0, 00). 

(3b) If < ^2 then < < = 'C2 '^nd neither b[ nor b'2 is continuous 
at Moreover, there exists mo G (— m*, m*) such that, as ^ I ^t, 

— mo mo + m* 

(2.20) 

(4) For 6 7^ &i(0'^2(0' the function m 1— > Q'^^{m) is minimized by a 
single number m — rri- {b, ^) G [— m*, m*] which satisfies 

--m^, ifb>bi{^), 
m-{b,0 { e (-rn*,m,), ^/^slO < b < 61(e), (2-21) 

- -m^, ifb<b2{0- 

(5) The function b 1— > m_(6, e) strictly increasing in b for b G 
[62(e)' continuous on the portion of the line b — bi{^) for 
which e > <^nd has jump discontinuities both along the line defined 
by b — 62(e) ^f^d along the portion of the line b = bi{^)for which < 
e < eu- There are two minimizers at the points where b ^ m_ (6, e ) 
is discontinuous with the exception of ib,^) = (6i(et),et) = (62(et);et) 
when et < eu. where there are three minimizers; namely, and mo 
from part (3b). 

As a simple consequence of the definitions, it is seen that the question of 
whether or not ei > e2 is equivalent to the question whether or not 

g{m^) - 2m^g'{m^) + ^-^(2m*)^5f"(m^) < g{-m^). (2.22) 

We claim that (2.22) will hold, or fail, depending on the values of the various 
parameters of the model. Indeed, writing e — tanh(ft;/2) we get 

g{m) = log(l + em) + const. (2.23) 

Regarding the quantity em as a "small parameter," we easily verify that the de- 
sired inequality holds to the lowest non-vanishing order. Thus, if m^ is small 
enough, then (2.22) holds for all k, while it is satisfied for all m^ whenever k 
is small enough. On the other hand, as k tends to infinity, 5'(m^) — g{—mi,) 
tends to log |^"'* , while the various relevant derivatives of q are bounded in- 
dependently of m+. Thus, as m* — > 1, which happens when J — > 00, the 
condition (2.22) is violated for k large enough. Evidently, the gap eu — et is 
strictly positive for some choices of J and k, and vanishes for others. 
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Since 61 (0) > 0, for ^ sufficiently small the ice region includes points 
with 6 > . Let us also show that the phase separation region can rise 
above 6 = 0; as indicated in the plot on the right of Fig. 1. Clearly, it suf- 
fices to consider 6 = and establish that for some J, k and ^, the absolute 
minimum of m t-* Qq,^{^) does not occur at ±m*. This will certainly hold if 

(Q-^)'(mO > and Ql^{-m^) > Qo,^{m^), (2-24) 
or, equivalently, if 

^ 2m ^ ^3'^^*^ ^^'^ - g{-m^)) > Wl (2.25) 

are both true. Some simple algebra shows that the last inequalities hold for 
some ^ once 

—^{g{m^) - gi-m^)) > 2m^g'{m^). (2.26) 

But, as we argued a moment ago, the difference g{m^) — g{~ni^^) can be made 
arbitrary large by taking k 1 and sufficiently close to one, while g'{m^) 
is bounded in these limits. So, indeed, the phase separation region pokes above 
the 6 = axis once k » 1 and J » 1. 

Comparing to the linear asymptotic of the phase diagram from Part I, we 
see that in the finite- volume system with minus (ice) boundary condition, the 
lines bounding the phase separation region are shifted upward and again are 
pinched together. In this case it is the line h = 61(C) that is parallel to its 
counterpart h — /i+(0)C for ^ > ^u, while h — 62(0 has the same asymptotic 
slope (in the limit ^ — > 00) as the function h — h'_{0)^. 



2.4. Properties of the spin marginal 

On the basis of Theorems 2.1-2.4, we can now provide a routine characteri- 
zation of the typical configurations in measure 

p±,ci,hi_ The following IS an 
analogue of Theorem 2.2 of Part I for the cases at hand: 

Theorem 2.4. Let d > 2 and let J > Jc{d) and k, > be fixed. 
Suppose that Assumption A holds and let (cl) and (hi) be two sequences such 
that cl > for all L and that the limits ^ and b in (2.5) exist and are finite. 
Let us define two sequences ofBorel probability measures pf^ on [— m^t, m*] by 
putting 

p±([-l,m]) = P^f '"^(Mi < mL"), m e [-1, 1]. (2.27) 

Then the spin marginal of the measure P^^^^ ''^ can again be written as a 
convex combination of the Ising measures with fixed magnetization; i.e., for 
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any set A of configurations {crx)xeAL' 

P^f ''^(^ X {0, 1}^-) = J ptidm) F^i{A\ML = L^L'^J). (2.28) 

Moreover, any (weak) subsequential limit of measures p\ is concentrated 
on the minimizers ofm ^ Q^^im). In particular, for b ^ ^2(0 limit 
p"*" = lim£,_>oo p\ exists and is simply the Dirac mass at m+ (6, ^ ) — the quantity 
from Theorem 2.2 — and similarly for p" — limi_>oo Pl and & 7^ ^2(0- 

On the basis of Theorems 2. 1-2.4, we can draw the following conclusions: 
For d-dimensional systems of scale L with the total amount of salt proportional 
to L'^~^ (i.e., the system boundary), phase separation occurs dramatically in 
the sense that all of a sudden a non-trivial fraction of the system melts/freezes 
(depending on the boundary condition). In hindsight, this is perhaps not so 
difficult to understand. While a perturbation of size L'^~^ cannot influence the 
bulk properties of the system with a single phase, here the underlying system 
is at phase coexistence. Thus the cost of a droplet is only of order so it is 
not unreasonable that a comparable amount of salt will cause dramatic effects. 

It is worth underscoring that the jump in the size of the macroscopic 
droplet at 6 = 61 or 6 = 62 decreases with increasing ^. Indeed, in the extreme 
limit, when the concentration is finite (nonzero) we know that no macroscopic 
droplet is present at the transition. But, presumably, by analogy with the re- 
sults of [4] (see also [3, 15]), there will be a mesoscopic droplet — of a particular 
scaling — appearing at the transition point. This suggests that a host of interme- 
diate mesoscopic scales may be exhibited depending on how cl and tend to 
zero with the ratio Iil/cl approximately fixed. These intermediate behaviors 
are currently being investigated. 

3. PROOFS OF MAIN RESULTS 

The goal of this section is to prove the results stated in Section 2. We begin 
by stating a generalized large deviation principle for both magnetization and 
the fraction of salt on the plus spins from which Theorem 2. 1 follows as an 
easy corollary. Theorem 2.2 is proved in Section 3.2; Theorems 2.3 and 2.4 are 
proved in Section 3.3. 

3.1. A generalized large-deviation principle 

We will proceed similarly as in the proof of Theorem 1.3.7 from Part I. Let 
A c Z'' be a finite set and let us reintroduce the quantity 
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which gives the total amount of salt on the plus spins in A. Recall that E^'"^ 
denotes the expectation with respect to the (usual) Ising measure with coupling 
constant J and plus/minus boundary conditions. First we generalize a couple 
of statements from Part I: 

Lemma 3.1. Let A C WJ^ be a finite set. Then for any fixed spin 
configuration a = (o-^) G { — 1, 1}"^, all salt configurations (Sx) G {0, 1}^ 
with the same A^a and Q\ have the same probability in the conditional mea- 
sure P^''''^{-\a = a). Moreover, for any S = {S^) G {0, 1}^ with Na = [c\A\\ 
and for any m G [—1, 1], 

P^'^'''{S occurs, Ma = L'^IAIJ) 

where the normalization constant is given by 

E 1{Ar,(5')=Lc|A|j}E^'-^(e'^'3-('^'^')+''^-('^)). (3.3) 
s'e{o,i}A 



Proof. This is identical to Lemma 1.3.3 from Part I. | 

Next we will sharpen the estimate from Part I concerning the total entropy 
carried by the salt. Similarly to the object ^^'^(cr) from Part I, for each spin 
configuration a = (a^) G {—1,1}^ and numbers ^, c G [0, 1], we introduce 
the set 

A'^^a) = {(S,) G {0,1}^: = lc\A\\, Ql = L^c|A|J}. (3.4) 

Clearly, the size of A^^^ia) is the same for all a with a given value of the 

magnetization; we will thus let A^^'^{m) denote the common value of \ A^/^{a) \ 
for those a with M^ia) = [m\A\\ . Let ^(p) = plogp + (1 — p) log(l — p) 
and let us recall the definition of the entropy function 

H(™..;c) = -l±^^(-^)-i^^(H(i^); (3.5) 

2 VI + m/ 2 \ 1 — m / 

cf formula (L2.7) from Part L Then we have: 

Lemma 3.2. For each r] > there exist constants Ci < oo and Lq < 

oo such that for all finite K d "L^ with |A| > Lq, all 9,c G [0, 1] and all m 
with \m\ < 1 — T) satisfying 

2^c ^ ^ , 2(1 -e)c ^ 

<l-ri and — < I - r] (3.6) 



1 + m 1 — m 
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we have 



IM 



< C,'^. (3.7) 



Proof. The same calculations that were used in the proof of Lemma 1.3.4 
from Part I give us 

with the substitutions Ma = L"^|A|J and Qa = L^c|A|J. By (3.6) and \m\ < 
1 — T], both combinatorial numbers are well defined once |A| is sufficiently 
large (this defines Lq). Thus, we can invoke the Stirling approximation and, 
eventually, we see that the right-hand side of (3.8) equals exp{|A|S(m, 6; c)} 
times factors which grow or decay at most like a power of | A| . Taking logs and 
dividing by |A|, this yields (3.7). | 

Our final preliminary lemma is concerned with the magnetizations out- 
side [— m^, m^] which are (formally) not covered by Assumption A. Recall the 
sequence of Wulff shapes Wl defined at the end of Section 1.2. Note that Wl 
contains, to within boundary corrections, L'^ sites. 

Lemma 3.3. Suppose that J > and let cl and Hl be such that Lcl 
and LHl have finite limits as L oo. For each e > 0, we have 

lim J-^logP±f ■'^^(iMiyJ > (m, + e)L^) = -oo. (3.9) 

L— >oo iv" i ^ ^ 



Proof. This is a simple consequence of the fact that, in the unadorned 
Ising magnet, the probability in (3.9) is exponentially small in volume — cf 
Theorem 1.3.1 — and that with LHl and Lcl bounded, there will be at most 
a surface-order correction. A formal proof proceeds as follows: We write 

P^f '"^(Ql = [ecLL% Ml = L^L'^J) = ^^I^pfl, (3.10) 

where 

KL{m, 9) = A'^^{m) e'^U-L'^J+«Lec,L''j pJ;^^(Mi = L^l'^J) (s.u) 

and where Yl is the sum of KL{m', 6') over all relevant values of m' and 6' . 
Under the assumption that both Iil and cl behave like 0{L^^), the prefac- 
tors of the Ising probability can be bounded between e~*-^^ and e"^^ , for 
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some C < oo, uniformly in 9 and m. This yields 

^_/3.12) 

The same argument shows us that Yl can be bounded below by e '^^ times 
the probability that My/^ is near zero in the Ising measure Inlightof J > 
Jc, Assumption A then gives 

lim inf — log > — oo. (3.13) 

On the other hand, by Theorem 1.3.1 (and the remark that follows it) we have 
that 

jim logP^^^^dM^^J > (m, + e)L<^) = -oo. (3.14) 

Plugging this into (3.12), the desired claim follows. | 

We will use the above lemmas to state and prove a generalization of The- 
orem 2.1. 

Theorem 3.4. Let d > 2 and let J > Jc{d) and k > be fixed. 
Let cl G [0, 1] and hi E be two sequences such that the limits ^ and b in 
(2.5) exist and are finite. For each m G [^m^,. m^] and 9 G (—1, 1), let Bl,^ = 
BlA^i, Cl, 9) be the set of all (a, S) G 1}^^'^ x {0, 1}^^ for which the 
bounds 

\Mw^ - mL'^l < eL'^ and \Qw^ - 9clL'^\ < eL'^'^ (3.15) 
hold. Then 

lim lim \ , ^' = -^bdm,9) + inf ^bdm',9'), (3.16) 

O'e[o,i] 

where ^^^(^{rn, 9) is as in (2.7). 

Proof. We again begin with the representation (3.10-3.11) for the 
choices HlL'^ ~ hL'^~^ and clL'^ ~ For m G [— m*, m^,] the last prob- 

ability in (3.1 1) can be expressed from Assumption A and so the only thing to 
be done is the extraction of the exponential rate of (m) to within errors 
of order o{L'^^^). This will be achieved Lemma 3.2, but before doing that, let 
us express the leading order behavior of the quantity S(m, 9; cl). Noting the 
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expansion y{p) — plogp — p + 0{p'^) for p J, we easily convince ourselves 
that 

2ecL 



E{m, 6; cl) = -Ocl flog - l) 

\ 1 + m / 



_(i_,K(,„gH(l^_i)+o(cD 

= cl-cl log Cl + CLT{m, 9) + 0(c|), 

where T(m, 6) is as in (2.4). (The quantity 0(c|) is bounded by a constant 
times c\ uniformly in m satisfying |m| < 1 —r] and (3.6).) Invoking Lemma 3.2 
and the facts that \ Wl\ — L'^ — 0{L'^~^) and Lc\ ^ as L ^ oo we now 
easily derive that 

A^^-(m) =exp{ri + L'^-^^T(m,^) + o(L'^-^)}, (3.18) 

where tl = —L\Wl\cl log(c/,/e) is a quantity independent of m and 6. 
Putting the above estimates together, we conclude that 

Xi(m, 9) = exp| tl - L^-^^b,^{m, 9) + o(L'^-i)} (3.19) 

where o{L'^~^) is small — relative to L'^~^ — uniformly in m G [— m^,mj,] 
and 9 E [0, 1]. It remains to use this expansion to produce the leading order 
asymptotics of Pp^'^''^'^{BL,e)- Here we write the latter quantity as a ratio, 

Pw^'^'iBL,) = ^^^^^pl^, (3.20) 

where KL^^(m, 9) is the sum of KL^m', 9') over all relevant values of (m', 9') 
that can contribute to the event Bl^^, while, we remind the reader, is the sum 
of KLimI , 9') over all relevant (m', ^')'s regardless of their worth. 

It is intuitively clear that the -factors in the numerator and denominator 
cancel out and one is left only with terms of order L'^~^, but to prove this we 
will have to invoke a (standard) compactness argument. We first note that for 
each (5 > and each (m, 9) G [— m^, m^] x [0, 1], there exists an e > and 
an Lq < oo — both possibly depending on m, 9 and 6 — such that, for L > Lq, 



\\og{KL,eim,9)e-'^) + ^f,,^im,9) 



< 5. (3.21) 



(Here we also used that ^b^^{m,9) is continuous in both variables 
on [— m*,m*] x [0, 1].) By compactness of [— m*,m*] x [0, 1], there exists a 
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finite set of (m^, 9k) 's such that the above e-neighboorhoods — for which (3.21) 
holds with the same S — cover the set [— m^, m^] x [0, 1]. In fact we cover the 
slightly larger set 

7^ = [-m, - e', + e'] x [0, 1], (3.22) 

where e' > 0. By choosing the e's sufficiently small, we can also ensure that for 
one of the k's, the quantity cS;,,^(mfc, 9k) is within S of its absolute minimum. 
Since everything is finite, all estimate are uniform in L > Lq on TZ. 

To estimate Yl we will split it into two parts, 1^ i and Yl^2, according 
to whether the corresponding (m', 9') belongs to TZ or not. By (3.21) and the 
choice of the above cover of TZ we have that log Yl,i is within, say, 36 of 
the minimum of (m, 9) i— > ^b^^{m, 9) once L is sufficiently large. (Here the 
additional S is used to control the number of terms in the cover of TZ.) On the 
other hand, Lemma 3.3 implies that Yl^2 is exponentially small relative to Yl^i. 
Hence we get 



lim sup 



-^log(r^e-'--)+ inf ^,,^{m',9') 



L 

6i'e[o,i] 



\m'\<_ 

Plugging these into (3.20) the claim follows by letting 5 [Q. | 



< 35. (3.23) 



Proof of Theorem 2. 1. This is a simple consequence of the compact- 
ness argument invoked in the last portion of the previous proof. | 



3.2. Proof of Theorem 2.2 

Here we will prove Theorem 2.2 which describes the phase diagram for the 
"liquid" boundary condition, see the plot on the left of Fig. 1. 

Proof of part (1). Our goal is to study the properties of the func- 
tion m I— (5^^(m). Throughout the proof we will keep J fixed (and larger 
than Jc) and write ^{■) instead of For m e [— m*, m*], let us define 

the quantity 

E^{m) = -^g{m) + J^{m). (3.24) 

Clearly, this is just Q^^{m) without the 6-dependent part, i.e., Q^^{m) — 
—bm + E^{m). Important for this proof will be the "zero-tilt" version of this 
function, 

E^{m) = E^{m) - E^{-m^) - (m + m,)D%^, (3.25) 
where D'^^ is the "slope of E^ between — and m*," see (2.11). Clearly, E^ 
and E^ have the same convexity/concavity properties but E^ always satisfies 

E^{-m^) = %{m^) = 0. 
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Geometrically, the minimization of Q^^{m) may now be viewed as fol- 
lows: Consider the set of points {{m, y) : y = E^{m)} — namely, the graph 
of E^{m) — and take the lowest vertical translate of the line y = bm which con- 
tacts this set. Clearly, the minimum of Q^^{m) is achieved at the value(s) of m 

where this contact occurs. The same of course holds for the graph y = E^{m) 
provided we shift b by D^^. Now the derivative E'^{m) is bounded below 
at m = — and above atm — rrii, (indeed, asm^ m^, the derivative diverges 
to — oo). It follows that there exist two values, — oo < bi{^) < &2(0 < oo> 
such that m = is the unique minimizer for 6 > bi{^), m — —rrii, is the 
unique minimizer for b < &2(0' neither m — 7X1^, nor m — —m^ is a 
minimizer when 62(0 < & < 

On the basis of the above geometrical considerations, the region where 61 
and 62 are the same is easily characterized: 

6i(^)=62(0 if and only if ^^(m) > e [-m^,m^]. (3.26) 

To express this condition in terms of ^, let us define T(m) = ^''{m)/g''{m) 
and note that E'^{m) > if and only if T(m) > Now, for some constant 
C = C( J) > 0, 

T{m) = C(m^ - m)-^ (m + cot(/^/2))^ (3.27) 

which implies that T is strictly increasing on [— m^t,m*) with T{m) — > 00 
as m t m^. It follows that either E'^ is concave throughout [— m*, m*], or there 
exists a T~^(^) e (— ?ti*, rrii,) such that i?^ is strictly convex on [—my,, T~^(^)) 
and strictly concave on (T~^(^),m*]. Therefore, by (3.26), < ^2(0 if 

and only if i?^(—m^) < 0, which is readily verified to be equivalent to { > 
This proves part (1) of the theorem. | 

Proof of parts (3) and (4). The following properties, valid for ^ > 
are readily verified on the basis of the above convexity/concavity picture: 

(a) For all 62(0 < ^ < there is a unique minimizer m^{b,^) of m 1-^ 
Q'l^{m) in [— m^,mj. Moreover, m+(6, ^) lies in (— m*,T~^(^)) and is 
strictly increasing in b. 

(b) For b = bi{^), the function m ^ Qb,^{'m) has exactly two minimizers, 
and a value mi (^) e (— m*, T~^(^)). 

(c) We have 62(0 = E'^{-m^)- 

(d) The non-trivial minimizer in (ii), nii (^), is the unique solution of 

E^{m) + (m* - m)E'^{m) = E^{m^). (3.28) 
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Moreover, we have 

6i(0=4(mi(0). (3.29) 

(e) As b tends to the boundaries of the interval (61 (0) ^2({)), the unique mini- 
mizer in (a) has the following limits 

lim m+(6, ^) = — and lim m+(6, ^) = mi(^), (3.30) 

where mi(,^) is as in (b). Both limits are uniform on compact subsets 
of (6,00). 

Now, part (3) of the theorem follows from (a) while the explicit formula (2.13) 
for 62(C) for C > Ct is readily derived from (c). For, ^ < ^t, the critical curve 
C I— > ^2(0 is given by the relation Q'^^{m^) — Q^^{—m^), which gives also 
the i < it part of (2.13). Continuity of 6 1— > m+(6, along the portion of 
& = &2 (0 foJ" > 'Ct is implied by (e), while the jump discontinuity at 6 = h\{^) 
is a consequence of (a) and (e). This proves part (4) of the theorem. | 

Proof of part (2). It remains to prove the continuity of &'i(0> identify 
the asymptotic of b[ as { — > 00 and establish the strict concavity of ^ i—> 61(C). 
First we will show that the non-trivial minimizer, mi(C), is strictly increasing 
with i. Indeed, we write (3.28) as F^{m) = 0, where F^{m) = E^(mi,) — 
E^{m) — {rrii, — m)E'^{m). Now, 

d 

—F^{m) = g{m) - g{m^) + {m^ - m)g'{m), (3.31) 
which is positive for all m e [— m*, m^) by strict concavity of g. Similarly, 

-^F^irn) = -E'l{m){m, - m), (3.32) 

which at m = mi (C) is negative because mi lies in the convexity interval of E^ , 
i.e., mi(C) e {—m^, T~^(C)). From (d) and implicit differentiation we obtain 
that m'i(C) > for C > Ct- By (3.29) we then have 

6'i = -^^^^ — ^^-^ (3.33) 
fn-k — TTii 

which, invoking the strict concavity of g and the strict monotonicity of mi, 
implies that b[{i) > 0, i.e., bi is strictly convex on (Ct, 00). 

To show the remaining items of (2), it suffices to establish the limits 



lim mi (C) = — m^ and lim mi(C) = m*. 



(3.34) 
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Indeed, using the fonner limit in (3.33) we get that b[{^) —g'{mi,) as ^ — > 
oo while the latter limit and (c) above yield that b[{^) b'^i^t) as { | which 
in light of the fact that bi{^) = 62(0 ^ 6 implies the continuity of b'^. 
To prove the left limit in (3.34), we just note that, by (3.28), the slope of 
at m = mi(^) converges to zero as ^ | Invoking the convexity /concavity 
picture, there are two points on the graph of m 1— > E^^{m) where the slope is 
zero: m^^ and the absolute maximum of E^. The latter choice will never yield 
a minimizer of Q'^^ and so we must have mi(^) — > as claimed. The right 
limit in (3.34) follows from the positivity of the quantity in (3.31). Indeed, 
for each m E [— m^,m^) we have F^{m) > once ^ is sufficiently large. 
Hence, mi({) must converge to the endpoint as { — >^ 00. | 



3.3. Remaining proofs 

Here we will prove Theorem 2.3, which describes the phase diagrams for the 
"ice" boundary condition, and Theorem 2.4 which characterizes the spin-sector 
of the distributions P^/^^ ''^. 

For the duration of the proof of Theorem 2.3, we will use the functions 
E^ and E^ from (3.24-3.25) with ^ = replaced by ^ = The 

main difference caused by this change is that the function m 1— > E^{m) may 
now have more complicated convexity properties. Some level of control is 
nevertheless possible: 

Lemma 3.5. There are at most two points inside [— m^,, m*] where the 
second derivative of function m 1— > E^{m) changes its sign. 



Proof. Consider again the function T{m) — ^"{m)/g"{m) which 
characterizes E'^{m) > by T{'m) > ^. In the present cases, this function 
is given by the expression 

T(m) = '^''}'^} = C{m^ + m)-^ (m + cot(/«/2))' (3.35) 
g [m) 

where C = C{J) > is a constant. Clearly, T starts off at plus infinity at m = 
— and decreases for a while; the difference compared to the situation in 
Theorem 2.2 is that T now need not be monotone. Notwithstanding, taking the 
obvious extension of T to all m > — m*, there exists a value mj e (— ?7t.*, 00) 
such that T is decreasing for m < rriT while it is increasing for all m > mf. 
Now two possibilities have to be distinguished depending on whether mx falls 
in or out of the interval [— m*, m*): 
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(1) mx > m^, in which case the equation T(m) = ^ has at most one solu- 
tion for every ^ and m i— > E^{m) is strictly concave on [— my^, T~^{^)) 
and strictly convex on (T'~^(^), m*]. (The latter interval may be empty.) 

(2) m-Y < m*, in which case the equation T{m) — ^ has two solutions 
for ^ e (T(mT), T'(myr)]. Then m i— > E^{m) is strictly convex between 
these two solutions and concave otherwise. The values of ^ for which 
there is at most one solution to T(m) = ^ inside [— m*, m*] reduce to 
the cases in (1). (This includes ^ = T(mT).) 

We conclude that the type of convexity of m i— > E^{m) changes at most twice 
inside the interval [— m*, m^] , as we were to prove. | 

The proof will be based on studying a few cases depending on the order of 
the control parameters and from (2. 18). The significance of these numbers 
for the problem at hand will become clear in the following lemma: 

Lemma 3.6. The derivatives E'^{m^) and E'^{mi,) are strictly increas- 
ing functions of^. In particular, for and ^2 defined in (2.18), we have 

(1) E'^{m^) <OifC< 6 and E'^{m^) >OifC> 6- 

(2) E'l{m,) <Oif^< 6 and E'^{m,) >Oif^> 6- 

Proof. This follows by a straightforward calculation. | 

Now we are ready to prove the properties of the phase diagram for minus 
boundary conditions: 

Proof of Theorem 2.3. Throughout the proof, we will regard the graph 
of the function m ^ E^ [m] as evolving dynamically — the role of the "time" in 
this evolution will be taken by ^. We begin by noting that, in light of the strict 
concavity of function g from (2.10), the value E^{m) is strictly decreasing in ^ 
for all m G (— m^, m*). This allows us to define 

= inf > 0: E^{m) < for some m e (— m^, my^)}. (3.36) 

Now for ^ = we have E^{m) > for all m E {—m^,,m^) while for ^ > 
^1, the minimum of E^ over (— m*,m*) will be strictly negative. Hence, we 
have < < 6- 

We will also adhere to the geometric interpretation of finding the mimiz- 

ers of m I— s> Q~^^{m), cf proof of part (1) of Theorem 2.2. In particular, for 

each ,^ > we have two values hi and 62 with 62 < &i such that the extremes 
— and are the unique minimizers for 6 < 62 and b > bi, respectively. 
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while none of these two are minimizers when b2 < b < bi. Here we recall 
that bi is the minimal slope such that a straight line with this slope touches 
the graph of at and at some other point, but it never gets above it, and 
similarly 62 is the maximal slope of a line that touches the graph of at — 
and at some other point, but never gets above it. 

As a consequence of the above definitions, we may already conclude 
that (1) is true. (Indeed, for ^ < |t we have E^{m) > and so the two slopes bi 
and 62 rnust be the same. For ^ > there will be an m for which E^{m) < 
and so 61 7^ 62.) The rest of the proof proceeds by considering two cases de- 
pending on the order of ^1 and ,^2- We begin with the easier of the two, ^1 > ^2^ 

CASE ^1 > ^2- Here we claim that the situation is as in Theorem 2.2 and, in 
particular, = ^i- Indeed, consider a ^ > ,^2 and note that E'^irrii,) > 
by Lemma 3.6. Since E'^{m) is negative near m = —m^ and positive 
near m — m*, it changes its sign an odd number of times. In light of 
Lemma 3.5, only one such change will occur and so [— m*, m*] splits into an 
interval of strict concavity and strict convexity of m h- > E^{m). Now, if is 
not equal ^1, we may choose ^ between and ^1 so that i?^(m*) < 0. This 

implies that E^{m) > for all m < rrii, in the convexity region; in particular, 
at the dividing point between concave and convex behavior. But then a sim- 
ple convexity argument E^{m) > throughout the concavity region (except 
at —rrii,). Thus E^[m) > for all m e (— m*, m*) and so we have ^ < It 
follows that It = ^1- 

Invoking the convexity/concavity picture from the proof of Theorem 2.2 
quickly finishes the argument. Indeed, we immediately have (4) and, let- 
ting = Ct> also the corresponding portion of (5). It remains to establish 
the properties of bi and 62 — this will finish both (2) and (3a). To this end we 
note that 61 is determined by the slope of E^ at rUi,, i.e., for ^ > 

6i(0=4K). (3.37) 

This yields the second line in (2.19); the first line follows by taking the slope 
of E(^ between — and m^. As for 62 > here we note that an analogue of the 
argument leading to (3.33) yields 

'm) = - "('"■)-/-'"♦) , ^ > |„ (3.38) 

mi -|- 

where mi = rni{^) is the non-trivial minimizer at 6 — &2(0- this case 
the argument analogous to (3.31-3.32) gives m[{^) < 0. The desired limiting 
values (and continuity) of b'^ follow by noting that mi(^) m^ as ^ I 
and mi(^) — >• — m* as ^ — >• 00. 
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CASE ^1 < Our first item of business is to show that < ^i- Consider 
the situation when ^ = .^i and m = nii,. By Lemma 3.6 and continuity, the 
derivative E'^^im^,) vanishes, but, since we are assuming < ,^2^ the second 

derivative E'^^ {m^) has not "yet" vanished, so it is still negative. The upshot 
is that is a local maximum for m i— E^^{m). In particular, looking at m 
slightly less than m^, we must encounter negative values of E^^ and, eventually, 
a minimum of E^-^ in (— m*, m^). This implies that ft < 

Having shown that Ct < Ci < C2, we note that for ^ e (ft, ^2)^ the func- 
tion m I— > E^{m) changes from concave to convex to concave as m increases 
from — to m*, while for ^ > ^2, exactly one change of convexity type oc- 
curs. Indeed, E^ is always concave near — and, when ^ < ^2, it is also con- 
cave at m^. Now, since ^ > ^t^ its minimum occurs somewhere in (— m^, m+). 
This implies an interval of convexity. But, by Lemma 3.5, the convexity type 
can change only at most twice and so this is all that we can have. For the 
cases ^ > ^2 we just need to realize that E^ is now convex near m — rrii, and 
so only one change of convexity type can occur. A continuity argument shows 
that the borderline situation, = ^2, is just like ^ > ^2- 

The above shows that the cases ^ > ^2 are exactly as for ^1 > ^2 (or, 
for that matter. Theorem 2.2) while C < Ct is uninteresting by definition, so 
we can focus on ^ e [^1,^2)- Suppose first that ^ > and let denote 
the interval of strict convexity of E^. The geometrical minimization argument 
then shows that, at & = bi, there will be exactly two minimizers, and a 
value mi(,^) G while at 6 = 62, there will also be two minimizers, — and 
avaluem2(0 ^ For 61 <b < 62, there willbe a uniqueminimizerm_(6, 
which varies between m2(C) and mi(^). Since E^ is strictly convex in I^, 
the map b 1— > m_(6, ^) is strictly increasing with limits mi(^) as 6 t ^i(0 
and m2(0 as 6 | ^2(0- Both nii and 7712 are inside (— m^, m*) so m_ under- 
goes a jump at both bi and 62- Clearly, mi(^) 7^ TO2(0 for all ^ e (^t, ^2)- 

At ,^ = ^t, there will be an "intermediate" minimizer, but now there is only 
one. Indeed, the limits of mi(,^) and m2{^) as J, {t must be the same because 
otherwise, by the fact that [mi(^), m2(0] is a subinterval of the convexity in- 
terval I^, the function E^^ would vanish in a whole interval of m's, which is 
impossible. Denoting the common limit by mo we thus have three minimizers 
at ,^ = ,^t; namely, ±m^ and rriQ. This proves part (4) and, letting = also 
part (5) of the theorem. As for the remaining parts, the strict concavity of bi 
and the limits (2.20) are again consequences of formulas of the type (3.33) and 
(3.37-3.38) and of the monotonicity properties of m-i and m2. The details are 
as for the previous cases, so we will omit them. | 

Proof of Theorem 2.4. As in Part I, the representation (2.28) is a 



26 



Alexander, Biskup and Chayes 



simple consequence of the absence of salt-salt interaction as formulated in 
Lemma 3.1. The fact that any subsequential (weak) limit of has all 
of its mass concentrated on the minimizers of is a consequence of The- 
orem 2.1 and the fact that m can only take 0{L) number of distinct values. 
Moreover, if the minimizer is unique, which for the plus boundary conditions 
happens when h ^ &2(0' subsequential limit is the Dirac mass at 
the unique minimum (which is m+(6, ^) for the plus boundary conditions and 
m_ (6, ^) for the minus boundary conditions). | 

ACKNOWLEDGMENTS 

The research of K.S.A. was supported by the NSF under the grants DMS- 
0103790 and DMS-0405915. The research of M.B. and L.C. was supported 
by the NSF grant DMS-0306167. 

REFERENCES 

1. K. Alexander, M. Biskup and L. Chayes, Colligative properties of solutions: I. Fixed 
concentrations, previous paper. 

2. K. Alexander, J.T. Chayes and L. Chayes, The Wulff construction and asymptotics of 
the finite cluster distribution for two-dimensional Bernoulli percolation, Commun. Math. 
Phys. 131(1990) 1-51. 

3. M. Biskup, L. Chayes and R. Kotecky, On the formation/dissolution of equilibrium 
droplets, Europhys. Lett. 60 (2002), no. 1, 21-27. 

4. M. Biskup, L. Chayes and R. Kotecky, Critical region for droplet formation in the two- 
dimensional Ising model, Commun. Math. Phys. 242 (2003), no. 1-2, 137-183. 

5. T. Bodineau, The Wulff construction in three andmore dimensions, Commun. Math. Phys. 
207 (1999) 197-229. 

6. T. Bodineau, D. loffe and Y. Velenik, Rigorous probabilistic analysis of equilibrium crys- 
tal shapes, J. Math. Phys. 41 (2000) 1033-1098. 

7. R. Cerf and A. Pisztora, On the Wulff crystal in the Ising model, Ann. Probab. 28 (2000) 
947-1017. 

8. R.L. Dobrushin, R. Kotecky and S.B. Shlosman, Wulff construction. A global shape from 
local interaction, Amer. Math. Soc, Providence, RI, 1992. 

9. R.L. Dobrushin and S.B. Shlosman, In: Probability contributions to statistical mechan- 
ics, pp. 91-219, Amer. Math. Soc, Providence, RI, 1994. 

10. P.E. Greenwood and J. Sun, Equivalences of the large deviation principle for Gibbs mea- 
sures and critical balance in the Ising model, J. Statist. Phys. 86 (1997), no. 1-2, 149-164. 

11. P.E. Greenwood and J. Sun, On criticality for competing influences of boundary and 
external field in the Ising model, J. Statist. Phys. 92 (1998), no. 1-2, 35^5. 

12. D. loffe and R.H. Schonmann, Dobrushin-Kotecky-Shlosman theorem up to the critical 
temperature, Commun. Math. Phys. 199 (1998) 117-167. 

13. R. Kotecky and I. Medved', Finite-size scaling for the 2D Ising model with minus bound- 
ary conditions, J. Statist. Phys. 104 (2001), no. 5-6, 905-943. 



Colligative properties of solutions 



27 



14. J.L. Lebowitz and C.-E. Pfister, Surface tension and phase coexistence, Phys. Rev. Lett. 
46(1981), no. 15, 1031-1033. 

15. T. Neuhaus and J.S. Hager, 2d crystal shapes, droplet condensation and supercritical 
slowing down in simulations of first order phase transitions, J. Statist. Phys. 113 (2003), 
47-83. 

16. R.H. Schonmann and S.B. Shlosman, Complete analyticity for 2D Ising completed, Com- 
mun. Math. Phys. 170 (1995), no. 2, 453-482. 

17. R.H. Schonmann and S.B. Shlosman, Constrained variational problem with applications 
to the Ising model, J. Statist. Phys. 83 (1996), no. 5-6, 867-905. 



